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Abstract 

An impact of integration over the paths of the Levy nights on the 
quantum mechanical kernel has been studied. Analytical expression 
for a free particle kernel has been obtained in terms of the Fox H- 
function. A new equation for the kernel of a partical in the box has 
been found. New general results include the well known quantum 
formulae for a free particle kernel and particle in box kernel. 
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1 Introduction 

In past decade it has been realized that the understanding of complex quan- 
tum and classical physic phenomena has required the implementation of the 
Levy flights random process pQ, [2|. It is well known that a position of 
diffusive particle increases as square root in time, x(t) ~ t 1 / 2 . The square 
root law is attribute of the Brownian motion model for diffusion. However 
for complex quantum and classical physic phenomena this temporal diffu- 
sive behavior has not been observed. Instead, more general evolution law 
x(t) ~ t 1 ' 01 with < a < 2 has been found. Thus, the well known diffusion 
law x(t) ~ t 1 / 2 is included as a special case at a = 2. The mathematical 
model to describe 1/a diffusion scaling is known as Levy nights. The Levy 
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flights is general framework to study 'unusual diffusion' in terms of fractional 
Fokker-Planck equation P, 

The scaling I /a law assigned to the Levy flights has been empirically 
observed in laser cooling of atoms [3], in ion dynamics in optical lattice [I], in 
anomalous transport [5] , in the measurement of the momentum of cold cesium 
atoms in a periodically pulsed standing wave of light [H] . The Levy flights are 
widely used to model a variety of physical phenomena such as kinetics and 
transport in classical systems, anomalous diffusion, chaotic dynamics, plasma 
physics, dynamics of economic indexes, biology and physiology, social science 
(see for example, [I], and references there). 

So too, the Feynman path integral ought to be generalized to describe 
complex dynamic phenomena at the quantum level. The Feynman path in- 
tegral is in fact the integration over Brownian-like quantum mechanical paths 
[7j, jH]. To generalize the Feynman path integral approach the integration 
has been expanded from Brownian-like to Levy flights paths [§], [TU] . 

In this paper we study an impact of integration over the paths of the 
Levy flights on the structure of a free particle quantum mechanical kernel. 
The kernel has been expressed in terms of Fox's if-function. 

The paper is organized as follows. We introduce path integral over the 
Levy paths in Sec. 2. It has been shown in Sec. 2 that the Levy path integral 
goes over into the Feynman path integral when self-similarity (or the Levy) 
index a=2. Thus, our new results go over into the well known quantum 
equations at the special case a = 2. 

Sec. 3 deals with a free particle quantum kernel. We have found Fox's 
if-function representation for a free particle quantum kernel. The represen- 
tation gives us an option to use many well established equations, identities, 
integrals involving the if-function for practical Levy path integral calcula- 
tions. 

Sec. 4 summarizes the Laplace and the energy-time transforms and the 
momentum representation of the Levy quantum mechanical kernel. 

New equation for the quantum kernel of a particle in a box has been 
found and discussed in Sec. 5. 
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2 Path Integral 



2.1 Feynman path integral 

If a particle at an initial time t a starts from the point x a and goes to a final 
point Xb at time £&, we will say simply that the particle goes from a to b and its 
trajectory (path) x{t) will have the property that x{t a ) = x a and x{t b ) = x b . 
In quantum mechanics, then, we will have an quantum-mechanical kernel, 
which we may write as K(x b t b \x a t a ), to get from the point a to the point b. 
This will be the sum over all of the trajectories that go between that end 
points and of a contribution from each. Following Feynman [7] we write the 
kernel in the form 

x{t b )=x b 

K(x b t b \x a t a ) = J Dx(r)exp{^S(x)}, (1) 

where S(x) is classical mechanical action considered as the functional of 
a particle trajectory x{t) 



S(r)= fdr[^--V(x(r),r))), (2) 



b / .2 

m x 



ta 



x(t b )=x b 

here V(x, t) is the potential energy and / Dx(t)... is defined as fol- 

£C (^a )— %a 

lows 



(3) 

Feynman's original path integral approach to non-relativistic quantum 
mechanics is based on the fundamental equations (JH)-©- 
Then the equation 



ip(x b ,t b ) = / dx a K(xbtb\x a t a )ip(x a ,t a ). (4) 
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gives the wave function ip(xb, tb) at a time t& in terms of the wave function 
if)(x a ,t a ) at a time t a . This fundamental equation describes the evolution of 
the quantum mechanical system in the terms of the wave function. 

2.1.1 The free particle 

To calculate the quantum-mechanical amplitude K^(xbtb\x a t a ) for a free 
particle, V(x) = we will follow Feynman [7j , 



x(t b )=x b ( _ ( 



b ■ 2 



K^{x b tb\x a t a )= J Dx{T)e^\ l -J dr^-Y (5) 

x{t a )=X a V *a 

With the help of Eqs.flU-© we have 



K^tM = *» {^) / n (^) " 2 x (6) 



-OO 



Jl eXp {We {Xj - Xj - l)2 }- 



The calculation is carried out as follows. It is easy to see that 1 



-1/2 



oo 



( 2iTihe\ f , ( im r, s2 , s2 ll 

I — I y cixiexp |— [(x 2 - xi) +(xi-x a )Jj = 

— oo 

Next we multiply this result by 

-1/2 



2niTie 



rn 



1 Here we use the definite integral 



dxe- ax ' 2+bx = J-e b2 / 4a 



4 



and integrate again, this time over x 2 - The result of integration is similar 
to that of Eq.(J7|), except that (x 2 — x a ) 2 becomes (x 3 — x a ) 2 and the expression 
2e is replaced by 3e in two places. Thus, we get 

In this way a recursion procedure is established which after iV — 1 steps 
gives 

(2mhNeY 1/2 f im . , 2 1 
(^^J eXP l^ (Xa ~ Xo) J" 

Since Ne = tb — t a , it is easy to see that the final result after N — 1 steps 

is 

^(0)/ + i + n ( 2iii?i(Jb-QY 1/2 l im(x b - x Q ) 2 | 

= ( — ™ — j - exp i 2^-o )• (8) 

The Eq.((HJ) implies 

Aa; oc f-V /2 (At) 1/2 , 

which means that space displacement Ax = Xb — x a and the time At = 
tb — t a scale are governed by the Brownian motion law. 

In terms of the Fourier integral the Feynman kernel K^\xbtb\x a t a ) given 
by Eq.(jHJ) can be written as 

K^\x b tb\x a t a ) = ^Jd P -e W - i At »~ h ta) ] ■ (9) 

— CO ^ ' 

2.2 The path integral over the Levy paths 

The "natural" way to generalize Eq.Q is [HI] 

vWt + \ + \ 1 1a (.P( x b-Xa) D a \p\ a {t b -t a )\ 
Kl J (x b t b \x a t a ) = — J dp -expi i i >, (10) 
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where D a is "fractional quantum diffusion coefficient" physical dimension 
of which is 

[D a ] = erg 1 "" • cm° • sec"", 

and a is the Levy index, and we consider 1 < a < 2. 
It is easy to see from Eq. (JlU)) that the scaling relation between space 
displacement Ax = x b — x a and the time increment At = t b — t a has a form 

l/a 



Ax oc (ft 0-1 D a ) (At) 1/a . 



This scaling relation implies that space displacement and the time incre- 
ment are governed by the Levy flights law. 

Using Eq.(fTU|) we can define the path integral in the phase space repre- 
sentation by the following way 



x{t b )=x b 



J Dx(t) J Dp(r)... = (11) 



oo 1 °° 

lim / dx 1 ...dx N _ 1 / dp^.ApNX 

N^oo J (27m) JV J 



Pi(xi-x a ) .D a e\p 1 \ a \ f .p N (x b -xat_i) D a e\p 



N\ 



a 



ex P |z jx...xexp|z % % J .., 

here again e = (t b — t a )/N. Then the kernel KL(xbtb\x a t a ) defined by 
Eq.((Tni) can be written as 

oo _^ oo 

K L (x b t b \x a t a ) = lim / dxi...dx N -i / dp x ...dp N x 



J i N I \ i N i N 

ex P | ^ _ J x ex P j-^^Xj \Pj\ a ~ ft £ Yl V ( x 3^ £ ) 

In the continuum limit iV — > oo, e —>■ we have 
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K L (x b t b \x a t a ) = 



(12) 



x{t b )=x b / _ t b 

J Dx(r) y Dp(r) exp j ~ y dr[p(r) x (r) - H a (p(r), x(t), t] 

x(t a )=X a y t a 

x(t b )=x b 

where the phase space path integral / Dx(r) JDp(r)... is given by 

2* (^a )— %a 

Eq.(|llj). x denotes the time derivative, H a is the Hamiltonian 

H a (p,x) = D a \p\ a + V(x,t) (13) 

with the replacement p — > p(r), x — > x(t) and {p(r),x(r)} is the particle 
trajectory in phase space. The hermiticity property of the Hamiltonian H a 
has been discussed in [TT] . 

The exponential in Eq.(|12j) can be written as exp{^S a (p, x)} if we intro- 
duce canonical action S a (p,x) for the trajectory pit), x(t) in phase space 

t b 

S a (p,x) = J dT(p(T) x (T)-H a (p{r)^{r),r)). (14) 

t a 

Since the coordinates xq, xn in the definition (fTTjl are fixed at their ini- 
tial and final points, xq = x a and xn = x b , the all possible trajectories in 
Eqs.llTTj). (|T2^1 satisfy the boundary condition x{t b ) = x b , x(t a ) = x a . We 
see that the definition given by Eq.(jllj) includes one more pj-integrals than 
Xj-integrals. Indeed, while x and x^ are held fixed and the x,- integrals are 
done for j = 1, iV — 1, each increment accompanied by one 

Pj-integral for j = 1, N. The above observed asymmetry is a consequence 
of the particular boundary condition. Namely, the end points are fixed in 
the position (coordinate) space. There exists the possibility of proceeding 
in a conjugate way keeping the initial p a and final p b momenta fixed. The 
associated kernel can be derived going through the same steps as before but 
working in the momentum representation (see, for example, [T2^). 

Taking into account Eq. (|lU|) it is easily to check on directly the consistency 
condition 

oo 

Kf\x h t h \x a t a ) = J dx'Kf\x b t h \x't') • Kf\x't'\Xata)- 

— oo 
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This is a special case of the general quantum-mechanical rule: the kernels 
for events occurring in succession in time satisfy 

oo 

K L (x b t b \x a t a ) = J dx'K L (x b t b \x't') ■ K L (x't'\x a ta), (15) 

— oo 

where KL(x b t b \x't') is defined by Eq. fT2|) . 

3 Fox H function representation for a free 
particle kernel Ki\xbh\x a t a ) 

Let us show how a free particle quantum mechanical kernel K^\x b t b \x a t a ) 
defined by Eq.(JTUJ) can be expressed in the terms of the Fox's if -function [T3] . 
[T4| . [To] . Apart from the quiet natural way in which the Fox's if -function 
enters in the path integral technique over the Levy flights, its derivatives and 
integrals are easily calculated by formally manipulating the parameters in 
the if -function. That is why the if-Fox's function is so important and useful 
for the path integration over the Levy paths. Note that if -function bears the 
name of their discoverer Fox ^3] although it has been known at least since 
1888 (see [Hj). 

To present a free particle kernel (x b t b \x a t a ) in terms of the Fox ff- 
function we accept the following strategy. Starting from the representation 
(fTU|) we obtain the Mellin transform of the quantum mechanical kernel. Com- 
paring of the inverse Mellin transform with the definition of the Fox's function 
|13j . [Tij yields the desired expression in terms of "known" function, i.e. Fox's 
if -function. 

Introducing for simplicity the notations 

X = x b - X a , T = t b - t a , (16) 

we rewrite Eq.(ll(J|) 

Ki } M = — J dp . exp U T -i— f —j. (17) 

— oo v ' 

One can see that the relation 

K i t\x^) = K { t\-x 1 r\ 
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holds. Hence, it is sufficient to consider K^\x,t) for x > only. 
Further, we will use the following definitions of the Mellin 



K (s, r) = J dxx s ~ x K L {x, r), (18) 
o 



and the inverse Mellin transforms 

c+ioo 



K l (x,t) = — f dsx- s K(s,r), (19) 

c—ico 

where the integration path is the straight line from c — ioo to c + ioo with 
< c < 1. 

The Mellin transform of the K^\x, r) is 

A °° 

Kf ] (s,t) = / dxx s ^K^\x,r) = 
o 

1 f , s-l f A j-P X ■ D a \p\ a T \ 

— - dxx dp - exp < i— i > . 

2irhJ J 1 ft ft J 

-oo v J 

By a change of variables 

' ft \ 1/a . / ft r 1A \ 



V — >■ £ — > ft 



JD a r J \iD a r J 

A 

, r) goes over to 

The integrals over d£ and cfc can be evaluated by using the equations 
(3.3) and (3.6) of the Ref. [IB]. Indeed, we have 

OO OO . , 

/ der 1 / ^exp{^-^} = -sin^— ^r( S )r(l - — ), (21) 

J J s — 1 2 a 
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where s — 1 < a < 2 and is the gamma function 2 . 
Inserting Eq. (J21|) into Eq. (J2(Jj) and using the functional relations for the 
gamma function, T(l — z) — — zT{— z) and Y(z)Y(l — z) = ix / sirniz, we find 



'(0) 



h V" 1 r( s )r(i^) 



Kl a\(h/iD aT y/«) r(±^)r(i±^ 

Then the inverse Mellin transform gives a free particle quantum mechan- 
ical kernel (x, r) 



c+ioo A 

K?\x,r) = ± J dax-K® 



J_I Th ( % V" 1 - r(s)r( ^ ) 

2m a' J aS [(h/iD a ry/<* % r(±f*)r(^)' 

c—ioo z z 

where the integration path is the straight line from c — ioo to c + ioo with 
< c < 1. Replacing s by — s we find 

KS\x,T)= (22) 

i / n \ - 1 j_ ~ c r ds (i(_n\ ^ \ s n-sn**) 

a\(h/iD a ry/ a ) 2m J \h\iD a r) ) r(±±*)rfe)" 

—c—ioo \ / z z 

The path of integration may be deformed into one running clockwise 
around R + — c. Comparison with the definition of the Fox's if-function (see 
Ens. (|?)5jl - (|Bn)l . Appendix A) yields 

K§\x,t)= (23) 

1 / H \ l/a x (1 - l/a, 1/a), (1/2, 1/2)" 
n\iD a r) X| (0,1), (1/2, 1/2) 

oo 

2 The gamma function T(s) has the familiar integral representation T(s) = J dtt s ~ 1 e~ t , 
Res > 0. 



a\(h/iD a ry/ a J ^ 
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Because of the Property 12.2.5 (see Appendix A) of the Fox's if -function 
the kernel K < //\x,t) becomes 



K f\x,r) = — E %2 



1.1 



ax 



l / h 



1/a 



h \ iD n r 



X 



(1, 1/a), (1,1/2) 
(1,1), (1,1/2) 



x>0. (24) 



Hence, for any x we have, 



Ki 0) (x,r) = —H, 



1 



ax 



i.i 

2.2 



1/a 



X 



h \iD a r J 

Remembering Eqs. (fTH|) . (|T7|) finally yields 

K^\x b t b \x a t a ) 



(1, 1/a), (1,1/2) 
(1,1), (1,1/2) 



(25) 



1 



a\x b - Xa 



1.1 
2,2 



h 



h \iD a (t b - t a ) 



1/0 \x -x I I (1 ' 1/2) ' 
16 (1,1), (1,1/2) 



(26) 



This is new equation for 1-D free particle quantum kernel (x b t b \x a t a ). 

Let us show that Eq. (}2T)j) includes as a particular case at a = 2 the well 
known Feynman quantum mechanical kernel given by Eq.fjHJ). Indeed, setting 
in Eq. (J26|) a = 2 and applying the series expansion Eq . (J63j) for the function 



#2,2 [| (^7) 



1/2 



I I (l,l/2),(l,l/2)i u 
X \ I (1,1), (1,1/2) J We naVe 



Kf\x,r)\ a=2 ' 



2^ V^2T, fc=0 

By substituting k —>■ 21 we obtain 



1/2 » / 1 / a \ i/2 \y 1 

M /' UvJ J fc! rc 1 -^- 



4°W)U 



1/2 



27i / 
Taking into account the identity 



Si r, \iD,r 



n \ 1/2 Y l x* 1 



(2i)\ra-i) 



(27) 
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^ , 1 \ , 1 . 7T 

+ *)r(- - 



2 2 COS7T2' 

and applying the Gauss multiplication formula 



we find that 



/24J-1 i 



(20!r(l-i) = ^| 7 (2) a, Z!. (28) 
(o). 



With the help of Eq.02EJ the kernel K\ ' (x, r)|o=2 can De rewritten as 



(Ph. -m 1 f n \ 1/2 ^ i ( n \ 1/2 \ 21 (-i) 1 * 21 



K L (x,r)| a=2 - 2 ^_ h ^— J g ^ - ^ 22//i 



(29) 



Since _D 2 = l/2m [5], JUj we finally obtain for the Feynman kernel 



^(».')-^.-)i- = S«p{w>- (30) 

It turns into Eq.(3-3) of Ref. [I] if we use Eq. (|TB|) . 

Thus, it is shown that a free particle Feynman kernel can be derived from 
the general equation (|2l)|). 

3.1 3-D generalization 

The above developments can be generalized to 3-D dimension. It is obviously 
that a free particle kernel given by Eq. ()10j) for the 3-D case has the form 
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Kf \r h t b \r a t a ) = J d 3 p ■ exp i i 



p(r 6 -r a ) .D a |p| a (t 6 - t a ) 



(2^)3y ^ i n & J ' 

(31) 

where r and p are the 3-D vectors. 

To express the kernel K^\r b t b \r a t a ) in the terms of the Fox's if -function 
we write 

Kf\r b t b \*ata) = 



OO 7T 2"7T s I . . . . 

1 2 fjafj \ .p\r b -r a \cosi> .D a \p\ a {t b - t a 



dpp 2 J d$ J dtp exp < z- 



(2nh) 3 J J J [ h h 

V ' -7T 



p|r b -r a | [ .-D Q |p| a (t 6 - t a ) 



f . pr 6 -r a I . 

y ctpp sm( ) exp < —i 



27r 2 ^ 2 |r 6 - r n | J h \ % 

With help of the formula 



. ,p|r 6 -r | 9 p|r 6 - r a | 

PSm( ^^ } = -^^-rj COS( ^^ } ' (32) 



(r b t b \r a t a ) can be written as 

Kf\v b t b \v a t a ) = -l-^-KP(x-,t b -t a )\ x=]r ^ Iab (33) 

where the kernel K^\x; t b — t a ) is 1-D kernel given by Eq. (|25jl . 

Let us note that Eq. ()33)l is just a special case of a general relation that 
holds between the D-mensional and D+2-mensional Fourier transforms of any 
isotropic function. An important consequence of Eq. (j33|) is that it allows us 
to evaluate the 3D kernel based on Eq. (|2"5j) for 1-D quantum kernel. Thus, the 
problem is to calculate the derivative of the 1-D kernel. Using the Property 
12.2.9 (see Appendix A) yields 

K ( °\r b t b \r a t a ) = (34) 
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1,2 

27ra |r 6 - rJ 3 ^ 3 ' 3 



1/a 



h \iD a (t b - t a ) 



rfe - r Q 



(1,1), (1, 1/a), (1,1/2) 
(1,1/2), (2,1) 



1,1 , 



This is new equation for a free particle quantum mechanical 3-D kernel. 
We see that in comparison with 1-D case the 3-D quantum kernel is expressed 

1 2 

in the terms of H 3 ' 3 Fox's if-function. 

Let's see that the 3-D kernel given by Eq. (j34j) goes over into the Feynman 
3-D kernel. Setting in Eq. (j34|) a = 2 and applying the series expansion 
Eq.fJOHJ) for the function H 33 we have 



Kf\r b t b \r a t a )\ a=2 



1 



E 



r(i + fc) (-i) fc 



n 



4nh\r b -r a \^ Q r{k)r(^) k\ \\iD a (t b -t a ) 
or 

(r b t b \r a t a )\ a=2 



1/2 \ j - A 

I [r 6 — r tt j 



(35) 



1/2 







E 



.i f_ 

h \i 



D 2 {t b -t a ) 



1/2 



r b - r a 



4vrft \iD 2 {t b -t a )) 8\r b -T a \^ k\Y(^) 

By substituting k — > 2/ and using Eq.(j2Hl) we observe that there exists 
the following relation 



7i 



iD2(t b -ta) 



1/2 



r 6 - r a 



l-fc 1 



which transforms Eq.()35jl into 



\ exp {-mD~(u-t a )}- (36) 



(v b t b \r a t a )\ a -_ 



Ti 



3/2 



4mhD 2 (t b -t a ) 



exp 



r& - rJ 



uhD 2 {t b - t a ) 



• (37) 
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Since D 2 = l/2m we got the Feynman 3-D quantum mechanical kernel 
K^°\r b t b \r a t a ) of a particle with the mass m (see, Problem 4-12, page 89 of 
Ref. 0), 

Thus, the general equation (J3~2j) includes the Feynman 3-D kernel as a 
special case at a — 2. 



4 Transforms of a free particle kernel 

4.1 The Laplace transform of a free particle kernel 

The Laplace in time transform k\ (x, s) of 1-D free particle kernel is defined 

as 

oo 

K { ?\x,s) = J dre- s -Kf\x,T), (38) 

Q 

where K^\x,r) is given by Eq. (J25j) . 

Using Eq.()38p and applying the series expansion for the function H 2 ' 2 
yield 



K?\x, s) = — £ n l^\ J-^ - (U^r) \x\ I / dre-r-^ 



1+fc 







(39) 



i ~ r(±±*)m - ±±*w-i) fe (i(hs\ 1,a , 



axs^ r(^)r(i^) fc! v ftV*A*, 

where we took into account that 

1 + fc. 



/i+fc , 
rfre- sr r-— = s" _1 r(l 



a 
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Applying the definition of the Fox's function (see Eqs. (f5HJ) - (f6T)j) . Ap- 
pendix A) we can write finally the Laplace transform K^\x, s) in terms of 
if -function 



l(_^V /a | I I (l,V«),(0,-l/a), (1,1/2) 
h\iD a ) (1,1), (1,1/2) 

(40) 

Putting in Eq.(|4U|) a = 2 and using the series expansion for the H 3 ' 2 - 
function we obtain the Laplace transform of a free particle kernel for the 
standard quantum mechanics 




4.2 The energy-time transformation 

In quantum mechanics an important role plays the Fourier transform of the 
kernel in the time variable, which is the fixed-energy kernel ki(x2,X\] E) 

oo 

k L (x 2 , xi; E)= J dheW^to-V ■ K^t^h), (42) 
h 

where Ki (a^ l^jii) is given by Eq. (fT2*j) . 

To make the integral convergent, we have to move the energy into the 
upper complex half-plane by an infinitesimal amount e. Then the Levy fixed- 
energy kernel becomes 

00 ih 

k L {x 2 , x 1 ;E) = J2 <Pn(x2)<P* n ( x i) ■ p f , • , ( 43 ) 

f^i E-E n + ie 

here 4> n (x) and E n are the eigenfunctions and eigenenergies of the Hamil- 
tonian H a defined by Eq.([13p. 

The small ie-shift in the energy E in (|43|) may be thought of as being 
attached to each of the energies E n which are thus placed by an infinitesimal 
piece below the real energy axis. When doing the Fourier integral ()42j) the 
exponential e -( l / h ) E ( t 2- t i) makes it always possible to close the integration 
contour along the energy axis by an infinite semicircle in the complex energy 



j3°W) = — ffg 

axs 
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plane, which lies in the upper half-plane for t 2 < ^1 an d in the lower half- 
plane for t 2 > h- The ze-shift guarantees that for t 2 < ti, there is no pole 
inside the closed contour making the kernel vanish. For t 2 > ii, on the other 
hand, poles in the lower half-plane give, via Cauchy's residue theorem, the 
spectral representation of the kernel ([43)1 . 

We see that the fixed-energy kernel k L (x 2 , x\\ E) and the kernel KL(x 2 t 2 \xibi) 
related each other by the inverse energy Fourier transform 



oo 

K^hlx^) = — ( dE e-WWto-V -k L {x 2 ,x^E). 
ItxTi J 



(44) 



Let us calculate a free particle fixed-energy kernel fc| , (r 2 , v^E) defined 
as follows 



oo 

fci 0) (r 2 ,r i; E) = J dt 2 e^ E ^ ■ Kf(r 2 t 2 \r ltl ). 



With help of Eq.flMJ) we have 



kf\r 2 , ri] E) 



1 



1 



2iia \r 2 — rJ 



(45) 



(46) 



H. 



1.2 

3,3 



1/a 



r 2 - ri 



(1,1), (1, 1/a), (1,1/2) 
h \iD a (t 2 - h) J rz ~ 1M (1,1), (1,1/2), (2,1) 

Using the Property 12.2.8 of if-function (see Appendix A) yields 
fcf(r 2 ,r i; £0 = — -J 1 



2iiaiE\r 2 — rj 



(47) 



rrl,3 
-"4,3 



£ \ 1/a 



h\ D r 



r 2 - r x 



(0, -1/a), (1,1), (1, 1/a), (1,1/2) 
(1,1), (1,1/2), (2,1) 



Thus the expression for the fixed-energy kernel k\ (x 2 ,T\]E) involves 
the iflf Fox's if-function. The inverse energy Fourier transform defined 
by Eq.(|44j) and the expression (|47|) allow to obtain the following alternative 
representation for the quantum mechanical kernel K^\r 2 t 2 \riti) 
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tfi 0) (r 2 * 2 |ri*i) 



n 



(27r) 2 m|r 2 — ri| 



-(i/h)E(t 2 -ti) 

dE x 

E 



d-1,3 
-"4,3 



1 / _EnV« _ (0, -l/«), (1, 1), (1, 1/a), (1, 1/2) 
n\ £>J |r2 111 (1,1), (1,1/2), (2,1) 



If we put a = 2 Eq. (J48j) goes over into the standard quantum mechanical 
fixed-energy kernel k^(r2,Yi,E)\ a= 2 = fc(°)(r 2 , i"i; E). Indeed, setting in 
Eq. (|48|) a = 2 and taking into account that in accordance with the definition 
of if -function (see Eqs.(|H5 |) -(|50 jl . Appendix A) 



H 



1,3 



1,3 



k ■ r 2 - ri 



(0,-1/2), (1,1), (1,1/2), (1,1/2) 
(1,1), (1,1/2), (2,1) 



- (k ■ |r 2 - ri|) 2 exp (-« ■ |r 2 - r x |) , 

we find the well known equation for the fixed-energy kernel (see, for ex- 
ample Eq.(1.390) at the space dimension D = 3, Chapter 1, [P2] ) 



A; (0) (r 2 ,r i;J E;) 



2m 



ATlTliK 2 \v2 — T\ 



m 



27r^i|r 2 — ri| 



- X (k • |r 2 - ri|) 2 exp (-/c • |r 2 - ri|) 
x exp (-k ■ |r 2 - n|) , 

1/q 



(49) 



where for simplicity the notation k = ^ ( — jy- J has been introduced. 
4.3 Momentum representation 

To find quantum kernel in the momentum representation let us introduce the 
momentum space wave function (p(p,t), 

<p{j>,t) = Jdre-^(v,t), 

where ip(r,t) is the wave function in coordinate representation 
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Then we consider the 3-D generalization of Eq. (0} 

i/)(r b ,t b ) = J dr a K L (r b t b \r a t a ) ■ ip(r a ,t a ). 

Substituting the equation for the wave function in coordinate represen- 
tation ip(r, t) in terms of the wave function in the momentum representation 
<p(p,t) yields 

f(Pb,h) = J dp a K L {p b t b \p a t a ) ■ (p{p a ,t a ), 

where the kernel in the momentum representation JCL(pbtb\Pat a ) is defined 
in terms of the kernel in coordinate representation KL(r b t b \r a t a ) as follows 

IC L (Pbt b \Pat a ) = J dr b dr a e-^ m+ ^ Ta ■ K L (r b t b \r a t a ). 

For example, for a free particle we have 

4°W&|P<A) = / dr b dr a e-r^n^ . K f\v b t b \v a t a ) = (50) 

= (2nh) 3 5(p a - p b ) ■ exp{--D a \p a \ a (t b - t a )}, for t b > t a . 

n 

We see from Eq. ()5()j) that a free particle kernel in the momentum space is 
expressed in terms of exponential function, while the kernel in the coordinate 
space has more complicated form. 

5 Particle in a box 

Now we are going to consider the impact of integration over the Levy flights 
paths on quantum kernel for a particle in 1-D box of length 2a confined by 
infinitely high walls at x = — a and x = a . From the eigenvalues [TU] 

E n = D a [—\ n a } Ka<2, (51) 



a 
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with the principal quantum number n = 1,2,3...., and corresponding 
eigenfunctions jTUj 

1 rara 

w^-^^Yi ( 52 ) 



it follows that the quantum kernel for a particle in the box has the form 

-^box(^FaO) = - V sin sin exp{--L> Q — ni}. (53) 

a a a n \ a J 

Here x a and x b are initial and final particle positions in the box. Note 
that the eigenfunctions ?p n (x) given by Eq. fjB^j) guarantee that the kernel 
satisfies the boundary conditions 

K hox (x b = a, t\x a , 0) = K haK (x b , t\x a = -a, 0) = 0, (54) 

enforced by the two infinite walls at x = —a and x = a at all times. 
Then Eq. (p)3*j) can be expressed as 

K hox (x b t\x a 0) = (55) 

l^f nn nix A i ( rrh\ a 

— > <^ cos — (x b - x a ) - cos — (x b + x a ) \ ex-p{--D a — n t\ 
2a ^ la a J n \ a J 

The next steps to transform the above equation are 

K hox (x b t\x a 0) = 



OO s 



f A xfP 7^ 7^ J rWfe-^Ji r 
2«ha , S 7 ^ - a 6XP[ ft ] - GXP[ T, 

L = — oo_ 



-oc 



exp{-^D a b| a 



1 ^ f , f r «p(^b - + 2/a) r -ip(a; 6 + a; a -2Za). 
2^ / aWexp[ ]-exp[ 



2yra , 



x 
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exp{--D a \p\ a t}. 

where the Poisson summation formula 3 has been applied. 
If we take into account the definition of a free particle kernel K^\xbt\x a 0) 
jiven by Eq.ffTU]) the kernel for the particle in the box Kh ox (xbt\x a 0) becomes 



oo 

K hm {x b t\xjd) = T. {K'°\x t + 2la,t\x a 0)-K 



l=—oo 



In terms of Fox's if -function K-b ox (xbt\x a 0) is 
K hox (x b t\x a 0) ' ^ 



a \{h/iD a ty/ a 



2la,t\x a 0)} . 



(56) 



(57) 



1,1 

2 



l=— oo 



1 



h \ iD n t 



I/a 



•l/a,l/a),(l/2,l/2) 
(0,1), (1/2, 1/2) 



#2,2 



1/q 



iD n t 



Xb — x a + 2la\ 



l/«,l/a),(l/2,l/2) 
(0,1), (1/2, 1/2) 



}• 



To understand Eq. fjBTj) for K hox (xbt\x a 0) from point of view of physics let 
us remind the analogy with the method of image charges in electrostatics. 
Following the method of image charges one can account the conducted plate 
by putting a negative charge complementing the positive charge. In this case, 
the image method yields an infinite number of charges of alternating signs. 
The original positive charge gives rise to two negative charges which are each 
an image corresponding to one of the two conducted plates. These images 
generate mirror images corresponding to other conducted plate and the image 
generation process has going on up to infinity. The quantum problem of 
particle in the box can be treated similar to the electrostatics problem of a 
charge between two parallel conducted plates. To apply electrostatic analogy 



a 

l— — oo 

where S is the Dirac delta function. 



l= — oo 
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to the quantum kernel let us start with simple physical situation of just 
one infinite wall and take a look at all paths running between x a and Xb 
in time t. As can be seen from the space-time diagram in Fig.l there is 
quantum path which does not cross the wall and which therefore contributes 
to the path integral. However between x a and Xb there exist also path which 
crosses the wall an even number of times. Since this path goes through the 
forbidden region, it should not contribute to the path integral. To eliminate 
contribution of those paths we have to apply a procedure to ensure that 
only quantum paths not crossing the wall will be taken into account. The 
procedure can be realized as follows. At first we write down a free particle 
kernel K^\xbt\x a 0) which disregards the wall. Then we have to subtract 
from K^\xbt\x a 0) the contribution of all paths which cross the wall. This 
can be done by constructing a free particle path integral with the same 
classical mechanical action as one which contributes to the Kj/ \xbt\x a 0). 
We take the original path up to the last crossing with the wall and then 
to continue along the mirror image of the original path. We thus end up 
at the mirror image -Xb of the original end point Xb- Note that a path 
running from x a to -Xb necessarily crosses the wall at least once. Therefore 
substracting a free particle kernel between x a and -Xb eliminates contribution 
from all paths which do not remain in the physical region x > 0. Finally we 
conclude that the kernel -ft'waii(^fe^l^aO) in presence of a wall can be expressed 
by substracting a free particle kernel going from x a to the mirror image -Xb 
from a free particle kernel going from x a to Xb, see Fig.l, 

K mll {x b t\x a 0) = K { °\x b t\x a 0) - K^\-x b t\x a 0). 

Expressing the kernel of a particle in the box in terms of free kernel works 
exactly in the same way. Indeed, a path intersecting both walls is subtracted 
twice, i.e. one time too often. Therefore, one contribution has to be restored 
which is done by adding another end point. Continuing the procedure one 
ends up with an infinite number of end points. The general rule to attribute 
a sign to each end point is that each reflection at a wall leads to factor -1. 
This procedure immediately gives us Eq. (j57|) . 

Using the way of consideration which leads from Eq. (fT7j) to Eq. fpIU]) we 
can see that at a = 2 the Levy quantum kernel Eq. (j57|) goed over into the 
well known Feynman kernel for a particle in a box (see, for example Eq.(6.19) 
and (6.20) in the |T2j). 
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6 Conclusion 



We have studied an impact of integration over the paths of the Levy flights 
on the structure of a free particle quantum mechanical kernel. Analytical 
expression of a free particle 1-D quantum kernel has been obtained in term 
of the Fox ^/-function. The 3-D generalization has been presented as well. 
The Laplace, energy-time and momentum transforms of a free particle kernel 
have been obtained and discussed. We have found the quantum kernel for a 
particle in the box. 

While the Feynman path integral is in fact integration over the Brownian- 
like paths the Levy path integral is the integral over the Levy flights trajec- 
tories. The path integral over the Levy flights generalizes at a < 2 and 
becomes at a = 2 the Feynman path integral. The new equations ()26p . ([34)1 . 
()40|) , PH|) , (|5U|) and (|57|) at the special case a = 2 go over into the well known 
quantum mechanical equations for a free particle kernel and a particle in the 
box kernel. 

7 Appendix A 
7.1 Fox i7-function 

In Sec. 3 we have expressed a free particle kernel in the term of Fox function 
H22 ■ Apart from the quiet natural way in which the Fox i7-function enters in 
the path integrals over the Levy flights, its fractional derivatives and integrals 
are easily calculated by formally manipulating the parameters in the H- 
function. That is why the Fox's if-function is important and useful for the 
Levy path integral calculations. 

Some properties of the if-function in connection with Mellin-Barnes inte- 
grals were investigated by Barnes [17] . Mellin jTH], Dixon and Ferrar [T^j. In 
an attempt to unify and extend the existing results on symmetrical Fourier 
kernels, Fox has defined the if-function in terms of a general Mellin- 
Barnes type integral. Asymptotic expansions and analytic continuations of 
the Fox function and its special cases were derived by Braaksma (201 ■ Many 
properties of the if-function are reported in the book ^3] along with appli- 
cations to statistics. 

Fox's if- function is defined by the Mellin-Barnes type integral [20] 
(we follow the notations of the book [Tl] ) 
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Tjm,n(^\ Tjm,n 



(dp, A p ) 
iKB q ) 



(58) 



Tjm^n 



(ai,Ai), (a p ,A p ) 
(b 1 ,B 1 ),...,(b q ,B q ) 

where function %(s) is given by 



h\ dszS x{s) 



and 



n r(bj - B jS ) n r(l - aj + A jS ) 
n r(i - bj + B jS ) n riaj-Ajs) 

j=m+l j=n+l 



z s = exp{sLog|2;| +iaigz}, 



(59) 



here m, n p and q are non negative integers satisfying < n < p, 1< m < q; 
and the empty products are interpreted as unity. The parameters Aj (j = 
l,...,p) and Bj (j = l,...,q) are positive numbers; aj (j = l,...,p) and bj 
(j = 1, ...,p) are complex numbers such that 



A j (b h + v)^B h (a j -\-l), 



(60) 



1, ...,n. 



for v, A = 0, 1, h — 1, m; j 

The L is a contour separating the points 



bj + v 
B~ 



, [j = l,...,m; z/ = 0,l,...), 



which are the poles of T(bj — Bjs) (j = 1, m), from the points 
/a, — v — l\ 



-4; 



(j = l,...,n; u = Q, 1,...), 



which are the poles of T(l — a^- + Ajs) (j = 1, n). The contor L exists 
on account of (|60p. These assumptions will be retained throughout. 

In the contracted form the H function in Eq. (J58|) will be denoted by one 
of the following notations: 
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H(z) 



TTm,n 



(a p , A p ) 
(b q ,B q ) 



The Fox if-function is an analytic function of z which makes sense (i) for 
every z ^ if /i > and (ii) for < \z\ < (3~ l if p, — 0, where 



/' X B J X - \i 

i=i j=i 



(61) 



and 



v q 



i=i j=i 



(62) 



Due to the occurrence of the factor z s in Eq. (j58|) . the if -function is in 
general multiple- valued, but is one- valued on the Riemann surface of log z. 

The H function is a generalization of Meijer's G function |21j . which is 
also defined by a Mellin-Barnes integral. The H function reduces to the G 
function if Aj = 1 and Bk = 1 for all j = 1, 2, ...,p and k = 1, 2, q, 



G™ q n (z) 



TTm,n 
PA 



{ax, 1), {dp, 1) 



1 (6 1 ,l),...,(6 ff ,l)J ' 
If further m — 1 and p < q, then the H function is expressible by 



H l,n 
p,q 



(ai,l)...(ap,l) 
(6i,l)...(6„l) 



EI T(l + 6i - aj) z b i 

3=1 



ftr(i + 6l _ 6 ) n r( % -6 1 ; 

j=2 j=n+l 



1 +6i - oi, 1 + 6i - a p ^ \ 
1 + 6 1 -6 2 ,...,1 + 6 1 -6, A ] )' 

in terms of generalized hypergeometric functions p F q p.4j. As far as many 
well-known special functions, such as error function, Bessel functions, Whit- 
taker functions, Jacobi polynomials, and elliptic functions are included in the 
class of generalized hypergeometric functions, all of them can be expressed 
in the term of the Fox's H function. 
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To represent of if -unction in computable form let us consider the case 

m . 

when the poles s = (bj + v)/Bj (j = 1, ...,m; v = 0, 1, ...) of T7 r(fej — 
are simple, that is, where 

B h (bj + \)±Bj(b h + v), j^h, 



h = l,...,m; is, X = 0, 1, 2, 

and the prime means the product without the factor j = h. Then we 
obtain the following expansion for the if -function 



m oo n r(i - aj + AjS hk ) n 'r(&j - Sjs hifc ) fc 



Kn n (z) = E E 



ft =i fc=0 n r(i-bj + B jShk ) n r( ai -A iSwfe ) fc - ^ 

j=m+l j=n+l 

(63) 

Sfc* = (6h + fc)/-B h , 

which exists for all z ^ if p, > and for < if /j = 0, where /j 

and (3 are given by Eqs. lJBTj) and (p32~j) . 

The formula (JH3|) can be used to calculate the special values of the if- 
function and to derive the asymptotic behavior for z — > 0. 



7.2 Some identities of the i-f-function 

The if -function possesses many interesting properties that are helpful for the 
calculations with the path integrals over Levy flights. We present the list of 
mainly used properties of the ii-function. The results of this section follow 
readily from the definition of the if -function given by Eq.(|58|) and hence no 
proofs are given here. 

Property 12.2.1 The if -function is symmetric in the pairs (a±, A±), (a n , A n ), 
likewise (a n+1 , A n+1 ), (a p , A p ); in (b u B ± ), (b n , B n ) and in (b m+1 , B rn+1 ), (b q , B q ). 

Property 12.2.2 If one of the (a,j,Aj) (j = 1, ...,n) is equal to one of 
the (bj, Bj) (j = m + 1, q) or one of the (bj, Bj) (j = 1, m) is equal to 
one of the (a,-, Aj) (j — n + 1, ■■■,p)], then the if-function reduces to one of 
the lower order, and p, q and n (or) m decrease by unity. 
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Thus, we have the following reduction formula: 



Tjm,n 



(ai, Ay), (op, Ap) 
(61, Si), (6,_i, (ai, Ai) 



(64) 



irm,n-l 



provided n > 1 and q > m. 
Property 12.2.3 



(a 2 ,A 2 ), (ap, A p ) 
(6i, J B 1 ),...,(6,_i, J B,_ 1 ) 



™ f I (ai,Ai),...,(a p ,A p ) 
™ ^ 1 (6 1 ,S 1 ),..,(6 9 ,S (I ) j 



(65) 



Tjn.m 



1 (l-6i, J B 1 ),...,(l-6„S,)\ 
z (1 - ai, Ai), (1 - a p , A p )y ' 



This is an important property of the if-function because it enables us to 

m n 

transform an if-function with \x = J2 Bj — J2 Aj > and arg x to one with 

i=i j'=i 

/i < and arg(l/x) and vice versa. 
Property 12.2.4 



\ H m,n ( I (ai,A 1 ),...,(a p ,A p )\ = ( k (a^kAJ, (a p ,kA p )\ 

k ™ [ 1 (6 1 ,S 1 ),..,(6 ff ,S ff ); ™ ^ 1 (6 1 ,fcS 1 ),..,(6 9 ,fcS 9 );' 

(66) 

where A; > 0. 
Property 12.2.5 

M<S::;S$)- 

/ 1 (a 1 + aA u A 1 ),...,(a p + aA p ,A p )\ 



p,q 



(b 1 + aB 1 ,B 1 ),...,(b q + aB q ,B q ) J' 



Property 12.2.6 
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n+i ( I (0, 7)5 •••) ( a pi A p ) \ m+lj „ / . (a p , Ap), (0, 7) 



TTm,n-ti I ( / ? ? \ jk ' P/ / 1 \r rrm-t-i,n 



(b q ,B q ),...,(r,>y)J v ' P+W V ' (rn),..,(b q ,B q ) 



where p < q. 
Property 12.2.7 



(68) 



rrm+l,n ( i ( a p> A>) , • • • > (1 _ r , l) \ _ / fiQ N 

p+w r' (i,7),..,(V5 ? ) J" (69) 



(1 - r,7), (a p ,A p ) 
(b q ,B q ),...,(l,j) _ 



where p < q. 

In the above Properties (12.1.2) - (12.1.6) the branches of the if-function 
are suitably chosen. 
Property 12.2.8 

d ^-v^-(^M<;;;^;:;t^)= (70) 

- aR m,n+l ( W_ | (1 ~ a,r), (oi,Ai), A p ) \ 

a P+1 ' 9 U r (6 1| S 1 ) I ... I (6 9 ,B 9 ) J' 

where m • n 7^ 0; r, p > 0; a, a* > 0; | argu;| < a*7r/2; Re a + r, min 

Reibj/Bj) > 0. 

Property 12.2.9 



^iCl" Ka^l (% '^ 



p,<? 



(71) 



a-v jjm 

As ± ± 1 



n+1 / ( n -y\$ I ( fl i> A i) 



where v is arbitrary, a, b > and a + (3 min(bj/Bj) > — 1 (1 < j ' < m). Here 
the qD v z notes fractional derivative of order v (see for definition Refs. 

123 )■ 
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8 Figure caption 

Fig.l. A path crossing the wall is cancelled by a path running to the mirror 
point of the end point 
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